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ABSTRACT 

The  high  Reynolds  number  (Re)  flow  past  a  rapidly  rotating  circular  cylinder  is  inves¬ 
tigated.  The  rotation  rate  of  the  cylinder  is  allowed  to  vary  (slightly)  along  the  axis  of  the 
cylinder,  thereby  provoking  three-dimensional  flow  disturbances,  which  are  shown  to  involve 
relatively  massive  (O(Re))  velocity  perturbations  to  the  flow  away  from  the  cylinder  surface. 

Additionally,  three  integral  conditions,  analogous  to  the  single  condition  determined  in 
two  dimensions  by  Batchelor  (1),  are  derived,  based  on  the  condition  of  periodicity  in  the 
azimuthal  direction. 
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1  Introduction 


One  of  the  most  important  and  fundamental  theorems  of  Fluid  Mechanics  to  have  been 
developed  in  relatively  recent  times  is  the  closed  streamline  theory  of  Batchelor  (1).  This 
states  that  in  a  two-dimensional  flow  involving  closed  streamlines 

y  curl  uj  ■  (Is  =  0.  ( 1 ) 

where  u  is  the  vorticity  and  ds  the  line  element  along  a  streamline  At  large  Reynolds 
numbers,  when  viscous  effects  can  be  neglected,  this  condition,  taken  with  the  inviscid 
(Euler)  equations  of  motion,  leads  to  the  result  that  the  vorticity  is  constant  within  a  region 
of  closed  streamlines. 

This  theorem  was  generalised  by  Blennerhassett  (2)  to  the  situation  which,  rather  than 
being  one  of  closed  streamlines,  involves  helical  streamlines.  Specifically  a  three-dimensional 
flow  with  three  velocity  components,  all  of  which  are  independent  of  the  axial  direction;  this 
implies  constant  axial  pressure  gradient.  This  analysis,  in  addition  to  confirming  Batchelor's 
(1)  result,  leads  to  an  additional  result  linking  the  constant  axial  pressure  gradient  and  the 
viscous  terms.  This  result,  taken  in  conjunction  with  the  inviscid  flow  equations  yields  the 
result  that  the  axial  velocity  must  be  proportional  to  the  streamfunction  for  the  motion  in 
the  plane  normal  to  the  axial  vector. 

Little  progress  appears  to  have  been  made  in  further  extending  the  ideas  of  Batchelor 
(1)  into  the  three-dimensional  regime,  in  spite  of  the  importance  of  flows  of  this  type.  It  is 
the  aim  of  this  paper  to  address  this  issue. 

Two-dimensional  flows  involving  closed  streamlines  involve  additional  interesting  sub¬ 
tleties.  Even  though  the  work  of  Batchelor  (1)  involves  the  inclusion  of  (small)  viscous 
effects,  there  is  insufficient  flow  physics  in  the  constant  vorticity  result  to  actually  determine 
its  value.  This  must  be  determined  by  recourse  to  boundary-layer  regions  located  on  body 
surfaces.  Riley  (3)  presented  an  example  of  two-dimensional  flow  inside  an  elliptic  container, 
driven  by  slippage  of  the  container  walls.  By  Batchelor’s  (1)  theorem  the  flow  in  the  core  of 
the  container  must  be  that  of  uniform  vorticity,  and  Riley  (3)  showed  how  just  one  particular 
value  of  the  vorticity  produced  the  appropriate  behaviour  (exponential  decay  at  the  outer 
edge)  inside  the  wall  boundary  layers.  These  ideas  have  subsequently  been  extended  by  other 
authors  to  other  situations  (see  below),  including  situations  involving  two  wall  layers,  see 
Duck  (4). 

One  of  the  classical  results  of  potential-flow  theory  is  that  of  uniform  flow  past  a  circular 
cylinder,  with  superimposed  circulation.  Physically  this  circulation  may  be  caused  by  the 
rotation  of  the  cylinder.  A  number  of  authors  have  studied  the  relationship  between  the  rate 
of  rotation  and  this  circulation.  Glauert,  (5)  considered  the  large  Reynolds  number,  large 
rotation  rate  problem,  whilst  Moore  (6)  considered  the  finite  Reynolds  number,  large  rotation 
rate  problem.  Lc~  (7)  and  Ingham  (8)  have  considered  fully  numerical  solutions  (finite 
Reynolds  numbers,  finite  rotation  rates),  whilst  Nikolayev  (9),  Negoda  and  Sychev  (10). 
Sychev  (11),  Lam  (12)  h  ave  investigated  the  large  Reynolds  number  (finite  rotation  rate) 
problem,  partly  using  the  ideas  of  Riley  (3).  There  appears  t,o  be  a  monotonic  relationship 
between  rotation  rate  and  circulation,  and  the  simple  model  of  outer,  potential  flow  together 
with  a  unidirectional  boundary  layer  <«  only  for  r»ta<ir»i  .-m...  dwr  <  rb  i<  a! 


value.  As  this  critical  value  is  approached,  a  stagnation  point  forms  off  the  surface  of  the 
cylinder  wall,  hut  inside  the  houndary  layer;  this  is  sufficient  to  disrupt  the  entire  model  at 
lower  values  of  rotation  rate. 

In  this  paper  we  consider  a  three-dimensional  analogue  of  the  above  problem.  We  take* 
a  uniform,  straight  circular  cylinder,  and  a  uniform  flow  far  Irom  the  cylinder  directed 
perpendicular  to  the  axis  of  the  cylinder,  and  the  cylinder  surface  is  rotating  with  a  laige 
angular  velocity,  which  is  dependent  on  axial  location.  We  further  suppose  that  the  Reynolds 
number  (as  defined  in  Sec  tion  2)  is  large'.  The  restriction  on  large  angular  velocity  renders 
the  problem  tractable  to  analytic  techniques,  whilst  pointing  the  way  in  whic  h  more-  general 
rotation  rates  may  be  tackled.  This  analysis  is  presented  in  Section  T  and  the-  corresponding 
numerical  results  are  in  Section  4.  In  Section  5  we  formulate  three  integral  conditions, 
analogous  to  the  condition  of  Batchelor  (1)  described  above.  In  Section  (>  we  present  our 
conclusions. 
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2  Formulation 


We  take  a  straight  circular  cylinder  of  radius  a .  together  with  polar  coordinates  (ur.O.dz). 
with  r  measured  radially.  ft  azimuthally  in  an  anti-clockwise  direction  and  c  axially.  The 
surface  of  the  cylinder  is  rotating  with  angular  velocity  ^*  -f  cos  Ac.  when*  *>  and  A  are 
constants.  Referred  to  the  coordinate  system  described  above,  we  suppose  that  far  from  the 
cylinder  the  fluid  velocity  takes  the  form  f  (cos  0,  —  sin  0.0).  The  flow  Reynolds  number  is 
defined  by 


v 

where  v  is  the  kinematic  viscosity  of  the  fluid  (assumed  constant).  In  this  paper  we  are 
primarily  interested  in  the  regime  Re  ^>1. 

Two  further  non-dimensional  parameters  may  be  defined,  namely 


n  = 


ua 

f’oo 


(3) 


and 


(4) 


describing  the  rate  of  rotation  and  degree  of  three-dimensionality  introduced  into  the  prob¬ 
lem. 

To  make  further  progress  in  this  paper,  we  assume  that  jr|  -C  1  and  further  that  jfl|  >  1 
(but  that  e  =  o(fl~} )).  Thus  we  focus  our  attent  ion  on  small  amounts  of  three-dimensionality 
and  high  rotation  rates. 

We  write  the  velocity  vector  as  f/^u,  and  the  pressure  as  where  p  is  the  density 

of  the  fluid  (assumed  constant).  We  then  develop  the  solution  in  powers  of  c,  namely 

u  =  u0  +  cu  +  0(f2),  (n) 

P  —  po  +  (p  +  0(t 2).  (b) 


For  r  —  1  =  0(1),  when  Re  3>  1  we  assume  the  basic  flow  takes  the  form 

u0  =  ^cos  ft  ^1  -  ,  -  sin  0^1  +  +  ^-,0^  ,  (7) 

this  being  merely  the  potential  flow  solution.  There  also  exists  a  boundary  layer  wherein 
r  —  l  —  0( Re-1/2).  The  problem  for  Uq  in  the  boundary  layer  has  been  studied  in  the  past 
by  a  number  of  authors,  as  noted  in  the  previous  section,  however  as  Q  —>  oc, 

u0  -»  (O.fi.O)  +  6>(!)  (8) 


for  r  —  1  =  0{  Re 

The  full  problem  for  u  is 


V  •  «  =  0, 

U  •  VU0  +  U„  •  Vfl  : 


T’  n  -L 


V?  2  - 


:? 


(9) 


(10) 


with  the  following  boundary  conditions 


(11) 

(12) 


u  =  (0,cosAr,0)  on  r—  i, 

u  — >  0  as  r  -t  oc. 

In  the  following  section  we  consider  the  solution  of  (10)  as  0  and  Re  — >  oo. 
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3  The  limits  Re  — >  oo,  Q  — >  oo 

Since  the  perturbation  velocity  u  is  triggered  by  an  0(1)  amount  (see  (11))  it  is  tempting 
to  speculate  that  this  in  turn  will  lead  to  general  perturbation  quantities  away  from  the 
cylinder  of  0(1).  However  this  is  certainly  not  the  case.  It  turns  out  that  the  only  plausible 
solution  for  r  =  0(1),  i.e.  the  only  asymptotic  development  that  leads  to  self-consistency 
involves  magnitudes  very  much  larger  than  0(1). 

Considering  first  the  solution  for  r  —  1  =  0(1),  then  on  account  of  the  linearity  of  the 
system  (9),  (10).  we  may  write 

u  =  (v(r,  9)  cos  Ac,  u(r,  9)  cos  \z,  u’(?\  9)  sin  Ac). 
p  =  p(r,8)cos\z. 

We  also  have  (see  (7))  that  as  Re  — »  oo, 

Uq  ( Vqi , fif/oo  +  Oqi, 0), 


(13) 

(14) 


where 


The  only  meaningful  series  development  of  the  perturbation  velocities  turns  out  to  be 


v  =  n  Re  u0(r)  +  ••••+  Re  V\  (r,  0)  -fl  •  •  • , 
u  =  Re  ui(r,  0)  +  •  ■  •  -fl  u2(r)  +  •  •  • , 
w  =  ft  Re  wo(r)  -f  •  •  •  -fi  Re  wi(t\ 0)  -fi  •  •  • , 

p  =  fl  Re  po (r.  0)  -fi  •  •  •  -fl  flpx  (v)  -fl  ■  •  • .  (17) 


These  expansions  prove  to  be  the  key  to  the  solution  of  our  problem,  although  a  prion,  it 
is  difficult  to  justify  these  expansions,  and  so  a  posteriori  verification  is  therefore  necessary. 
In  Section  5,  some  prima  facie  justification  for  this  form  of  solution  development  is  given. 

Consider  now  the  leading  order  terms  in  (17)  that  are  independent  of  0.  It,  turns  out, 
perhaps  surprisingly,  that  the  flow  is  not  governed  by  the  inviscid  Euler  equations,  but 
rather  is  predominantly  viscous  in  nature.  Taking  0(H)  terms  in  the  radial  momentum 
equation  yields  the  following  equation 


_  dpx 
r2  dr 


+  Vorr  H - I'or  —  A  Vq 


Vo 

~r2  ' 


The  0(Re  ’)  terms  in  the  azimuthal  momentum  equation  yield 


1 

U2rr  H - u2r 

r 


A a«2  -  '4  =  0, 


(18) 


(19) 


whilst  the  0{£l)  terms  in  the  axial  momentum  equation  yield 


»’Orr  4-  -U’0r  -  A2l«o  =  -  (20) 

/• 

The  continuity  equation,  to  leading  order,  may  be  written 

vqt  +  ~i’o  +  Au)0  =  0.  (21 ) 

r 

With  this  system  there  is  no  difficulty  in  imposing  the  full  and  correct  boundary  condi¬ 
tions;  we  therefore  have 


Vo(r  =1)  =  u-'0(r  =  1)  =  0, 
u2(r  =  I)  =  1, 


with 


wQ,  wq,  u-2  —>  0  as  ?•  — >  co. 

The  solution  for  u 2  can  be  written  in  terms  of  Bessel  functions,  namely 

tf.(Ar) 

U‘~  A'i(A)  ’ 

whilst  it  is  possible  to  eliminate  u>0  and  pj,  from  (18),  (20),  (21)  to  yield 

2  m  /  3 


(22) 

(23) 

(24) 


tttt  2  tti  /  i  _  %  It 

v0  +  ~v0  -  +  2A  J  yo 

/  3  2A2Vy 


+ 


2 A2  3  \  2u.)A2 


,.3 


v0  +  I  A  +  7T  -  v° 


(25) 


to  which  the  following  additional  (alternative)  boundary  condition  is  appropriate 

v'0(r  =  1)  =  0.  (26) 

Once  u2  is  determined,  it  is  quite  straightforward  (although  a  numerical  task)  to  deter¬ 
mine  Vo(r),  Wo(r);  numerical  results  will  be  presented  in  the  following  section.  It  is  worth 
noting  that  ?>o,  u-2  and  wo,  as  described  above,  also  represent  an  exact  solution  of  the  lin¬ 
earised  Navier  Stokes  equations  ((9),  (10))  in  the  absence  of  the  uniform  flow. 

We  now  turn  to  consider  the  leading  order  terms  in  (17)  that  are  dependent  upon  both 
r  and  0  and  which  are,  indeed,  determined  primarily  through  invisrid  equations.  Again,  the 
linear  nature  of  (9),  (10)  greatly  simplifies  the  solution  technique,  this  time  allowing  us  to 
write 


t>i  =  Mt(r)sin0, 

U]  —  «i(r)ros0, 
a > i  =  Till  (?•)  sin  0. 

Po  ~  pu(r)  cos  0.  (27) 
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We  may  also  write 


V'oi  =  V'oi  cos  0, 

(for  —  ^ot  sin  9, 

(28) 

with 

i 

H 

3 

1^ 

<'oi  -  — ( 1  H — j)- 

vi 

(29) 

Taking  terms  0(Q  Re)  in  the  radial,  azimuthal  and  axial  momentum  equations  rexper- 

tively,  leads  to 

—  dv  o  dV  oi  v,  2ui  (tf>0 

*  01  dr  +V°  dr  +  H  ”  r’  ~  ~  dT’ 

(20) 

dVoi  1  _  ,  f;oi  1 

v0  A  +  V0  ~  ~Po 

dr  r 1  r  r 

(31) 

T7  dw0  1 

V°>  dr  +  -  XP^ 

(32) 

whilst  continuity  leads  to 

Id.  ,  ill 

-  ~(rvt) - -  +  Au?i  -  0. 

r  dr  r 

(33) 

On  this  system  we  may  impose  decay  of  all  components  of  the  solution 
v  —  1,  we  may  have  only 

as  r  — >  oo,  whilst  on 

v,(r  =  1)  =  0. 

(34) 

The  no-slip  constraints  on  ¥|  and  ui[  are  therefore  violated,  but  these 

may  be  rectified  by 

the  inclusion  of  a  thin  boundary  layer  of  thickness  0(Sl~'i  Re-*)  on  r  ~ 

=  1. 

Specifically,  we  write 

Y  -  (r  -  1)^2  Re 2  =  0(1), 

«i  -  Ui(Y)eiB  +  e.c^ 

(35) 

(30) 

1! 

+ 

<*«v 

(37) 

and  so 

in,  =  u | y v  -  /p0(  1 ), 

(38) 

mc,  =  w\yy  +  Ap0(  1 ), 

(39) 

giving 

«i  =  -p„(l)[l  -c-(,+l,V/^], 

(10) 

ihj  =  — ?Ap„(l)  |l  _  f-C+OV/x/Ij 

(ID 

We  see  therefore,  that  on  the  r  —  1  —  0(1)  scale,  intriguingly  both 

viscous  and  iuviscid 

effects  are  important  simultaneously,  and  to  a  large  degree,  independent  tv.  In  the  following 
section  we  consider  a  number  of  numerical  results  arising  from  the  results  of  tins  section, 
and  go  on  to  consider  the  limits  A  — >  oo  and  A  — >  0. 


4  Results  and  large/small  A  behaviour 

The  system  (22),  (25),  (30)-(34)  was  solved  using  a  conventional  fourth  order  Runge-Kutta 
method.  Results  are  shown  in  figure  !  for  ]i0(r  =!)(---  -).  j>\  (r  —  1)  (solidus).  UjT(r  = 
— ),  where  the  dependence  of  these  quantities  with  A  is  shown. 

We  may  make  some  further  analytic  progress  by  considering  the  limits  of  large  and  small 
A.  Taking  first  A  — ►  oo,  then  by  (24) 

— f’*\  (42  j 

where 

V'  =  A(r  -  1)  =  0(1).  (43) 

i.e. 

«2r|r_,  — *  —A.  (44) 

The  above  implies  the  perturbation  to  the  flow  is  confined  to  within  a  thin  r  —  1  =  0(A~'  ) 
boundary  layer  (although  we  must  impose  the  restriction  that  A  =  o(fiO  Re^)  in  order  that 
this  layer  remains  outside  the  Y  =  0(1)  layer  discussed  in  t he  previous  section).  After  some 
algebra,  it  is  possible  to  show  that 

(45) 

(46) 

(47) 

(48) 

Comparison  of  (47)  evaluated  on  V  =  0  with  the  corresponding  numerical  results  of  figure  1 
shows  good  agreement  as  A  increases.  Unfortunately  it  is  not  possible  to  compare  pQ(r  —  1) 
with  />„( V  =  0)  without  further  substantial  algebra  since  p0(V  =  0)  =  0(A~4). 

The  alternative  limit  of  A  — ►  0  is  slightly  more  complicated,  because  two  key  radial 
lengt hscales  emerge.  For  r  —  0(1).  we  have  t hat  (24)  reduces  to 

u 2  ~  — .  (44) 

r 

wliirh  interestingly  is  the  azimuthal  velocity  component  corresponding  to  a  line  vortex,  i.e. 
effectively  the  two-dimensional  result.  I  here  is  also  a  large  radial  scale.  R  =  Ac  =  ()( 1 ).  (We 
do  have  the  restriction  that  A  =  o(il)  in  order  that  over  the  lengt  hscales  under  consideration, 
/’oo  remains  dominant  over  I  «i  in  (16)). 

For  r  —  0(1),  the  solution  develops  in  the  form 

t’o  =  AJ  [log  Af’o(r)  +  (){ 1 )]  .  (50) 


v  o  = 


w0  - 


j»Ve  +0{X  1 

Ll(iv^ir)+r;,A-’) 


V ' 


Pi 

Po 


A 

-Y 


3 


--V  +0(A-2). 


1  ,  1  V-3  :)l>2  0 


,  .  -r‘  +  —  y*  - 

A3  V8  12  8 


-Y 

8 


0(A"4; 


<s 


with  the  leading  order  general  solution  taking  the  form 

,  I) 

t’o  =  Ar  +  Hr  log  r  +  Cr  d - .  (  >1  ) 

r 

For  R  —  0(1),  we  must  consider  the  first  two  leading  terms  in  the  solution  development, 
namely 

i’o  =  A  [log  Ay0(/?)  +  i’\[R)  +  O((log  A)-1)]  .  (o2) 

The  complete  analytic  solution  for  r0(ff)  and  V\(R)  does  not  appear  possible,  however  it 
suffices  to  consider  the  limit  R  — >  0,  for  which 

v0{R)  ~  A0R  +  Bo/flog  R  +  C0R ;J  +  (fid) 

H 

whilst  MR)  (which  is  forced  directly  by  the  uj,  term  on  the  right-hand-side  of  (25)  )  takes 
the  form 


MR)  ~  A\ft+  fi, /? log «  +  r, ft3 

+  ~  -  jfl(logK)2. 

Boundedness  constraints  demand 


C  -  Do  =  Di  =  0, 


whilst  matching  of  (49)  as  r  — *  oo  with  (51)  as  R  — ►  0  rec|uires 


Bo  = 
B  = 


Further  imposing  the  two  boundary  conditions  on  r  —  1  yields 


vq 


-r - v  log  r  — 

8  4” 


8r 


This  then  leads  to 


Wo  =  -A  log  A  log  r  +  0(A), 
Pi  =  -1/r2. 


(54) 


(55) 


(56) 

(57) 


(58) 


(59) 

(60) 


Fortunately  it  is  possible  to  determine  the  leading-order  term  of  p i  without  substantial 
algebra  (and  additional  numerical  effort),  even  though  p\  is  an  order  lower  in  log  A  than  may 
first  be  expected;  equation  (4.19)  yields  the  result  that  p\{r  —  1 )  — >  —  1  as  A  — >  0,  a  result 
that  agrees  with  our  numerical  results.  However  it  is  not  possible  to  determine  p0(r  —  1) 
without  further  substantial  algebra  and  numerical  effort. 

In  the  following  section  vve  go  on  to  consider  results  for  more  genera!  classes  of  two- 
dimensional  flow’s  involving  closed  streamlines,  which  are  perturbed  in  some  three-dimensional 


manner. 


5  The  periodicity  requirement  and  the  associated  in¬ 
tegral  conditions 

The  resuit  of  Hatch  -lor  (11  may  he  viewed  as  arising  Item  a  condition  of  periodicity  within 
the  region  of  closed  streamlines  of  the  various  physical  flow  quantities.  It  is  needed  lieeauso 
the  Ruler  equations  are.  themselves,  not  sufficient  to  enforce  periodicity:  an  alternative 
viewpoint  is  that  tin  Fnler  equations  are  iinalde  to  capture  the  mean  flow  physics,  i.e.  the 
flow  corresponding  to  zero  wavenumber,  a  result  that  is  little  surprising.  1  his  point  is  equally 
important  in  the  three-dimensional  context,  and  is  now  investigated. 

We  may  write  (')).  ( 10)  in  the  form 

U0  A  Cj  -f  U  A  =  V  //  -f  He-1  V  A  T.  { 6  1  ) 

where  II  is  the  0(c)  total  lead.  i.e. 

H  -  p  +  u0  •  u.  (62  ) 

and  we  have  written  the  vorticity  vector 

a-’  —  u.’o  +  (J  +  Off2).  (6.1) 

Following  Batchelor  (1)  and  Blennerhassett  (2)  the  work  of  this  section  is  most  efficiently 
carried  out  in  terms  of  a  coordinate  system  based  on  the  two-dimensional  .  undisturbed 
st reamlunrt ion.  specifically  in  terms  of  ((.••£.  c).  Here  i/  represents  the  st reamfunction  of 
the  undisturbed  (ut»)  flow,  and  £  is  orthogonal  to  v  and  c:  c  remains  the  axial  coordinate. 
If  the  undisturbed  flow  Uo  is  irrotational.  then  £  may  be  taken  to  be  the  standard  velocity 

potential  o.  1  he  infinitesimal  line  element  is  ( — .h2d£.dz).  Here  h2  is  tlie  £  coordinate 

\  %  i 

metric,  which  in  the  case  of  an  irrotational  u0  is  merely  1 /<•  and  %  —  juoj. 

Referred  to  the  (i.\£.c)  coordinate  system,  the  velocity  vector  may  be  written 

u  =  ( c  v  cos  Ac.  (/,,  4-  <  u  cos  Ac.  tic  sin  Ac)  4-  0(c2 ),  (64  ) 

and  the  pressure  as 

P  —  Pa  4-  t pi  cos  Ac  +  ()((*).  (6o) 

Referred  to  the*  (i.\£.c)  coordinate  system,  we  may  write  the  vorticity  terms  as  follows 


—  (  0.  0.  o.’()  ) 

=  (0.0.  -p 

h  2  OV 


=  (u.’i  sin  Ac.wVi sin  Ac, cos  Ac) 

=  f  — j — ( tr<:  4-  An/;2).sin  Ac(  —  Ac 


c/o  cos  Ac 

h2 


() 

Jh 


(uli. 


(I 


<7o"V) 


(66) 


(67) 
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We  may  then  write  the  three  components  of  (til)  in  the  following  form 


IIuJq 

Vo 


+ 


~  t-’U-’o  — 


uhj^,  !  ()  I  r 

t>2  <K  \<h, 

[a*, 

Vo  '  l>2<lo  <K 

1  Op  1  0 

T~^7  +  7~ttt('/o») 


2  Vo 

<V/i  Re 

+ 


A/; 


( (is  i 


+  Re  1  |  A^q  -  (/()—“ 
Oir 


Ap  Re-1 

- 1 — - — 

%  hi 


'V  (Af  V'/o 


(69) 


ru  t 


The  continuity  equation,  in  terms  of  these  variables  is 

#  . ,  ,  0  f  u\  Xwh, 

5-r<**-)+jF(-)  + - =«- 


dip 


(K  Wo, 


Vo 


■i) 


if  we  integrate  each  of  (68) -(71 )  around  a  complete  circuit  in  £  (lying  entirely  within  a  region 
of  closed  streamlines),  then  periodicity  of  the  flow  demands 


hi 


UU>0  Ilk  2 

— j-  +  uhit  ~  Voe - 


Vo 


Vo 


hi  0£ 

Vo  Oir 


Re 


-1 


Vo 


~W 


—  A  I)  2^2 


<K  =  o, 


vu)q  -f-  Re  1 

-A/3  Re'1  0 
—  + 


0u)% 

Au-'i  —  r/o  "7 


h2<K  =  0. 


(hr 

{h-iui-i) )  h2(K  =  0. 


(72) 

(78) 

(71) 


Vo  hi  (hi 

These  results  are  quite  general,  and  indeed  exact  within  the  framework  of  the  linearised 
Navier  Stokes  equations  (9),  (10).  If  we  now'  direct  our  attention  on  the  problem  considered 
earlier  in  this  paper,  then  uo  is  irrotational.  and  so  in  line  with  our  comments  regarding  ( 
above,  we  can  write 


Op 


“J  |^V0v  M  +  7 ~  +  R< 


-1 


OijT\\  A 

~iK 


u-2 

Vo 


<K  =  n. 


/ 


AW) 


L  Vo 


Oir 


<l(  =  o. 


(77)) 


(76) 


X}>  '  ”  f)  (h 2^2))  K  =  «■ 


Vo  Oir 

Note  that  the  result-  of  Batchelor  (1)  is  retrieved  from  (76)  by  allowing  A  — >  0. 


At  this  stage  we  can  parth  linn  the  solution  development  of  Section  T  Rquation  (77) 
above  .strongly  suggests  that  u/2  =  0  (Re),  which  in  turn  suggests  that  e  and/or  w  must 
also  be  0(  R<  j.  This  is  partly  confirmed  by  (7.')).  Indeed,  it  can  be  shown  that  our  solut  ions 
given  ill  Section  3  do  satisfy  the  above  conditions,  which  clearly  illustrate  the  importance  of 
viscosity  on  the  r  —  1  —  0(1)  scale. 
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6  Conclusions 


In  this  paper  we  have  considered  the  efforts  of  t  liree  dimensionality  ini  rodneed  into  the 
problem  of  a  rapidly  rotating  circular  i  ylinder  in  a  uniform  llow.  The  rllecis  in<  hide  a. 
substantial  (0[  Re))  response  in  the  t  liree  veloeit  v  nnnpoiienls.  in  the  bulk  ol  the  fluid .  1  his 
effect  is  confirmed  by  t he  analysis  ol  Section  n  in  which  t liree  integral  conditions  were  d *  i  ived 
(analogous  to  the  solitary  integral  condition  in  Batchelor's  (1)  two-dimensional  work).  AH 
tli<’  indications  are  that  this  massive  response  with  t liree  dimensionality  will  be  a  generic 
feature  of  similar  flows.  It  is  likely  that  the  results  of  this  paper  have  repercussions  for 
important  practical  applications,  perhaps  the  most  important  of  which  being  that  of  high 
lift  ae  rofoils. 

This  paper  has  deliberately  focused  on  the  large  rotation  rate  problem  (i.e.  U  x  u 
However  the  extremely  important  regime  of  12  r=  0(1)  remains  to  be  studied.  1'his  regime 
involves  a  number  of  additional  questions,  perhaps  the  key  aspect  being  the  nature  of  the 
solution  as  the  closed  streamline  region  of  the  base  (two-dimensional)  flow  is  exited.  It  could 
well  be  that  some  form  of  mild  discontinuity  exists,  for  which  a  thin  shear  layer  would  be 
required.  This  whole  problem  would  be  a  non  trivial  numerical  undertaking,  but  it  is  to  he 
expected  that  the  O(Re)  velocity  scales  will  persist.  Kqually  we  have  implicitly  assumed 
here  that  12  ^  0(Re»).  n  positive  integer,  although  if  this  condition  is  relaxed,  it  seems 
likely  that  the  solution  will  he  modified,  but  in  a  relatively  minor  manner. 

The  <  =  ()( 1 )  problem,  at  this  stage,  would  appear  to  be  a  formidable  task,  quite  possibly 
involving  a  solution  of  the  full  Navier  Stokes  equations.  Finally,  although  our  study  has 
concentrated  on  sinusoidal  disturbances,  in  I  lie  axial  direct  ion.  it  is  of  course  st  raight  forward 
to  extend  our  ideas  to  non  periodic  axial  disturbances  using  standard  transform  techniques. 
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